Abstract In normal mesoscopic metals of a ring topology persistent currents can be 
induced by threading the center of the ring with a magnetic flux. This phenomenon 
is an example of the famous Aharonov-Bohm effect. In the paper we study the cur- 
rent vs the external constant magnetic flux characteristics of the system driven by 
both the classical and the quantum thermal fluctuations. The problem is formulated 
in terms of Langevin equations in classical and quantum Smoluchowski regimes. We 
analyze the impact of the quantum thermal fluctuations on the current-flux character- 
istics. We demonstrate that the current response can be changed from paramagnetic 
to diamagnetic when the quantum nature of the thermal fluctuations increases. 
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1 Introduction 

Mesoscopic systems lay at the border between macroscopic and microscopic worlds. 
Mesoscopic systems, consisting of a large number of atoms, are too big to study their 
properties by the quantum methods of individual atoms, and are too small to apply 
physical laws of the macro-world. For their description and modeling one should 
combine both methods appreciating and recognizing the role of quantum and classi- 
cal processes. The appearance of a new length scale - the phase coherence length l c of 
electronic wave functions - introduces various regimes for transport phenomena in- 
fluenced by quantum interference effects. The mesoscopic regime is characterized by 
small length scales and low temperatures. When the temperature is lowered, the phase 
coherence length increases and the mesoscopic regime is extended to larger length 
scales. At sub-Kelvin temperatures, the length scales are of the order of microme- 
ters. The most prominent mesoscopic effects are: the Aharonov-Bohm oscillations 
in the conductance of mesoscopic structures, the quantum Hall effects, the universal 
conductance fluctuations and persistent currents in mesoscopic normal metal rings 
threaded by a magnetic flux. Persistent currents have been known to exist in super- 
conductors in which they are related to the existence of a state with zero resistance 
and the fact that a superconductor is a perfect diamagnet. The existence of persistent 
currents in normal (i.e. non-superconducting) metals is, in fact, a manifestation of 
the famous Aharonov-Bohm effect. Persistent currents have been predicted by Hund 
in 1938 1 1 J and re-discovered later by others [2 |. Inspired by these findings, the first 
experiment was performed in the early 1990s by measuring the magnetization of an 
array of about ten million not connected micron-sized copper rings |3). Other experi- 
ments also supported existence of persistent currents [4 |. We should also recall recent 
definitive measurements of persistent currents in nanoscale gold and aluminium rings 
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Q6]. The team [5| has developed a new technique for detecting persistent currents 
that allows to measure the persistent current over a wide range of temperatures, ring 
sizes, and magnetic fields. They have used nanoscale cantilevers, an entirely novel 
approach to indirectly measure the current through changes in the magnetic force it 
produces as it flows through the ring. The second team (6) has studied thirty three 
individual rings, in which they have employed a scanning technique (a SQUID mi- 
croscope). The rings are very small, each only between one and two micrometers in 
diameter and 140 nanometers thick. They are made of high-purity gold. Each was 
scanned individually, unlike past experiments on persistent currents conducted by 
other groups. In total they were scanned approximately 10 million times. Both works 
mark the first time that the theory has been experimentally proven to a high degree. 

In our previous papers [7| we have proposed a two-fluid model for the dynamics 
of the magnetic flux that passes through a mesoscopic ring. It is described in terms 
of an ordinary differential equation with an additional random force. It is analogous 
to the well known model of a capacitively and a resistively shunted Josephson junc- 
tion (8.'| . The classical part consists of 'normal' electrons carrying dissipative current. 
The quantum part is formed by those electrons which maintain their phase coher- 
ence around the circumference of the ring (it is a counterpart of the Cooper pairs of 
the electrons in the superconducting systems). The effective dynamics is than deter- 
mined by a classical Langevin equation |9 1 with a Johnson noise describing classical 
thermal equilibrium fluctuations. For low temperatures, quantum nature of thermal 
fluctuations should be taken into account. To this aim we apply the approach based 
on the so called quantum Smoluchowski equation as introduced in Ref. [ 10 1 and in 
other versions in the following Refs. lllllll2l[T3l . 

The paper is organized as follows. First, in the Sec. 2, we present a model of ca- 
pacitively and resistively shunted Josephson junction in order to demonstrate readers 
the analogy between both models. Next, in the Sec. 3, we briefly present our model 
for the flux dynamics in the normal metal rings. In the Sec. 4, we define the quantum 
Smoluchowski regime following by the Sec. 5, presenting the dimensionless variables 
and parameters. In the Sec. 6, we study the current characteristics in the stationary 
states for both classical and quantum Smoluchowski domain. We end this work with 
the summary and conclusions. 



2 Superconducting ring 

For clarity of modeling the current characteristics in non-superconducting rings, we 
present the well-known approach to describe a quasi-classical regime of supercon- 
ducting rings. To this aim, let us consider a superconducting loop (ring, cylinder, 
torus) interrupted by a Josephson junction. This element is a basic unit of various 
SQUID devices. The phase difference ip of the Cooper pair wave function across the 
junction is related to the magnetic flux <j> threading the ring via the relation J5] 



ip = 2n(n - (p/(f> ), 



(1) 
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where 2nn is the phase change per cycle around the ring and <fio = h/2eis the flux 
quantum. When the external magnetic field is applied, the total flux is 

- ^ e + LI, (2) 

where (j> e is the flux generated by an applied external magnetic field, L is the self- 
inductance of the ring and / is the total current flowing in the ring. We model the 
Josephson element in terms of the resistively and capacitively shunted junction for 
which the current consists of three components [ 14 1, namely, 

I = I c + Ir + Ij=^^, (3) 

where Ic is a displacement current accompanied with the junction capacitance C, Ir 
is a normal (Ohmic) current characterized by the normal state resistance R and Ij is 
the Josephson supercurrent. In the right-hand side, the relation |2]) has beed used. 

Combining Eqs. ([TJ-([3j with the second Josephson relation dip/dt — 2eU/K, 
where U is the voltage drop across the junction, we get the Langevin-type equation 
in the form 1 8 1 



d 2 ^ .Idj . 0-^ e , /2fc B T 

c m> + RTt + Iosm(/) = —fr + V -ir m - (4) 

This equation has a mechanical interpretation: it can describe the 'position' <p(t) of 
the Brownian particle moving in the washboard potential 

W{4>) = ^ ^L^ 2 - J ° cos ^ (5) 

The first term originates from the external bias and the self-inductive interaction 
of the magnetic flux whereas the second term is the supercurrent modified by the 
quantum flux. The ubiquitous thermal equilibrium noise r(t) consists of Johnson 
noise associated with the resistance R. The parameter ks denotes the Boltzmann 
constant and T is temperature of the system. The Johnson noise is modeled by 6- 
correlated Gaussian white noise of zero mean, (r(i)) — 0, and unit intensity, i.e., 

(r(t)r(u)) = s(t-u). 



3 Normal metal ring 

Now, let us consider a non-superconducting ring. When the external magnetic field is 
applied, the actual flux is given by 

<t> = 4>e + LI, (6) 

where (f> e is the flux generated by an applied external magnetic field, L is the self- 
inductance of the ring and / is the total current flowing in the ring. At zero tem- 
perature, the ring displays persistent and non-dissipative currents Ip run by phase- 
coherent electrons. It is analogue of the Josephson supercurrent Ij, At non-zero tem- 
perature, a part of electrons becomes 'normal' (non-coherent) and the amplitude of 
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the persistent current decreases. Moreover, resistance of the ring and thermal fluc- 
tuations should be taken into account. Therefore for temperatures T > 0, the total 
current consists of three parts, namely, 

I = I C + Ir + Ip = ^^. (7) 

What we need is the expression for the persistent current I p. It is s function of the 
magnetic flux <f> and depends on the parity of the number of coherent electrons. Let 
p denotes the probability of an even number of coherent electrons and 1 — p is the 
probability of an odd number of coherent electrons. Then the persistent current can 
be expressed in the form ifTSIl 

Ip=I P (<l>)=pI E (ct>) + {l-p)Io{4>), (8) 

where 

oo 

I E {(f>) = Io(<f> + 0o/2) = Io A n (T/T*)cos(nk F l) sm(2mr<^ ), (9) 

n=l 

where Iq is the maximal current at zero temperature. The temperature dependent 
amplitudes are determined by the relation ITT31 

A (T/T*) 4T cxp(-nT/n 

A n (2 II ) - ^ 1 _ cM _ 2nT/T ^ , (10) 

where the characteristic temperature T* is proportional to the energy gap Ap at the 
Fermi surface, kp is the Fermi momentum and I is the circumference of the ring. If 
the number N of electrons is fixed then kp = irN/l and the persistent current takes 
the form 

oo 

I P (cj>) =I J2 MT/T*) sin(2mr0/0 o )[p + (-1)"(1 - p)}. (11) 

71=1 

As a result, from Eq. (j7]i we obtain the equation of motion in the form Q 



d 2 ct> 1 dcf> 1 2k B T 



where the "potential" V(tf>) reads 



V&) = ~ (0 - 4> e f + M £ A "f /r * } cos ( 2mr J ) [p + (-!)"(! 



2L vr ' ™ U ^-J 2n7r 

n=l 

In the above Langevin equation, C and L are, respectively, the capacitance and induc- 
tance of the ring. It was shown [ 16], that the energy associated with long-wavelength 
and low-energy charge fluctuations is determined by classical charging energies and 
therefore the ring behaves as it were a classical capacitor. The flux dependence of 
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these energies yields the contribution to the persistent current ifTTl . The capacitance 
becomes essential if the ring accommodates a stationary impurity or a quantum dot 
1 18 1. Moreover, in the mesoscopic domain the standard description of a capacitor in 
terms of the geometric capacitance (that relates the charge on the plate to the volt- 
age across the capacitor) gives way to a more complex notion of capacitance which 
depends on the properties of conductors 1 19 1. 

Note that Eqs. Q and ( 12 1 have a similar structure. The difference is not only in 
the form of the potential but also in temperature dependence of the potential in the 
case the normal metal rings. 



4 Quantum Smoluchowski regime 

In both models (for superconducting and non-superconducting rings), thermal equi- 
librium fluctuations are modeled as classical fluctuations of zero correlation time. 
When temperature is lowered, quantum nature of fluctuations starts to play a role, 
fluctuations become correlated and leading quantum corrections should be taken into 
account. It is not a simple task and a general method how to incorporate quantum 
corrections in a case described by Eq. ( fT~2| > is not known. However, in the quantum 
Smoluchowki regimes 1101 . where the charging effects (related to the capacitance 
C) can be neglected, the system can be described by the "overdamped" Langevin 
equation - the so-named quantum Smoluchowski equation 110111 11 . For a Brownian 
particle it corresponds to neglecting inertial effects related to the mass of a particle. 
The quantum Smoluchowski equation has the same structure as a classical Smolu- 
chowski equation, in which the diffusion coefficient Do = ksT / R is modified due 
to quantum effects like tunnelling, quantum reflections and purely quantum fluctua- 
tions. In terms of the Langevin equation ( p"2} , it assumes the form 



This equation has to be interpreted in the Ito sense 1201 . The modified diffusion co- 
efficient Dj\(tf>) takes the form IfTTl 

D *M = l-AV%)/k B T > (15) 

where the prime denotes differentiation with respect to the argument of the function. 
The quantum correction is characterized by the parameter A, It measures a deviation 
of the quantal flux fluctuations from its classical counterpart, namely, 

A = (0 2 ) Q - (0 2 ) c , (16) 

where (■ ■ •) denotes equilibrium average, the subscripts Q and C refer to quantal 
and classical cases, respectively. Let us determine the range of applicability of the 
quantum Smoluchowski regime. The classical Smoluchowski limit corresponds to 
the case when charging effects can be neglected. Formally, one can put C = in the 



7 



inertial term of Eq. (fT2j, which is related to the strong damping limit of the Brownian 
particle. In the case studied here it means that [10] 



uj CR < 1, 



(17) 



where the frequency ujq is a typical frequency of the bare system and its inverse 



corresponds to a characteristic time of the system. In such a case, Eq. ( 16 1 takes the 
form Q 



A 



fiR 



7 + ^1 + 



(18) 



2TrCRk B T / 

where the psi function &(z) is the logarithmic derivative of the Gamma function and 
7 ~ 0.5772 is the Euler constant. 

The separation of time scales, on which the flux relaxes and the conjugate observ- 
able (a charge) is already equilibrated, requires the second condition, namely, 



ujqCR <C ksT/huio- 



In the deep quantum regime, i.e. when 



27rCi?' 



the correction parameter ( 18 i simplifies to the form 



A : 



TiR 



7 + In 



2-KCRk B T 



(19) 



(20) 



(21) 



In order to identify precisely the quantum Smoluchowski regime, we have to deter- 
mine a typical frequency ujq or the corresponding characteristic time tq oc 1/ljq. 
There are many characteristic times in the system, which can be explicitly extracted 
from the evolution equation ( 12 1, e.g. CR, fr/ksT, 4>q/{RIq). The characteristic time 
r = L/R is the inductive time of the ring and for a typical mesoscopic ring, L/R 
is in the picosecond range. Therefore, in the quantum Smoluchowski regime, all the 
above inequalities (17i, (19i and (20i should be fulfilled for ujq — 2tt/tq. Because 



the diffusion coefficient cannot be negative, the parameter A should be chosen small 
enough to satisfy the condition Da{4>) > for all values of <j>. 

The "overdamped" Langevin equation ( p"4} describes a classical Markov stochas- 
tic process and its probability density P(4>, t) obeys the Fokker-Planck equation |20|, 
namely, 



1 d_ 
R dt 



P(cj>,t) = 



d 



dV(4>) 
d<t> 



P(0,t) 



S 2 
defy- 



[D A {4,)P[ct>,t)] 



(22) 



We wish to analyze an averaged stationary current (J) flowing in the ring which can 
be obtained from Eq. Q: 

(i) = \ m - m , (23) 

where the averaged stationary magnetic flux (0) is calculated from the equation 

(<f>)= </> PtfW, P(0)= lim P(fct), (24) 
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where P(4>) is a stationary probability density, 

P(<f>) = CoD^icj)) exp [-Um , (25) 

where Cq is the normalization constant and the generalized thermodynamic potential 
U(4>) reads 

U(<f>) = J dV ^D- A 1 ^) dcf>. (26) 

Because the potential V(<f) depends on the external flux <p e , the averaged stationary 
current ( |23] > is a non-linear function of <j> e , Eqs. ( |23] > - p6j ) form a closed set of 
equations from which the current characteristics (I) = f{<p e ) as a certain function / 
of the external magnetic flux <j> e can be obtained. It is an analogous of the current- 
voltage characteristics for electrical circuits. 



5 Dimensionless variables and parameters 

To analyze the current-flux characteristics in the stationary state, we first introduce 



dimensionless variables and parameters. The rescaled flux x — 4>/4>q. Then Eq. (23 1 
can be rewritten in the dimensionless form 

i = (x) - x e , i = (7)L/0o, x e = e /<t>o, (27) 

where i, (x) and x e are dimensionless averaged current, averaged flux and external 
flux, respectively. 

The stationary probability density p(x) takes the 

p(x) = N D- 1 (x)e X p[-* x (x)], (28) 

where No is the normalization constant and the generalized thermodynamic potential 
reads 

dV(x) 
dx 



!P A (aO= / a ^>D-\x)dx. (29) 



The rescaled potential reads 

1 



V(x) = -(x-x e ) 2 + B{x), (30) 



where 



00 A (T } 

B ^ = a E C0S ( 2mra )[P + - P)] W 

n—l 

with the dimensionless temperature To = T/T* and a = LIo/<f>Q. The rescaled 
modified diffusion function D\(x) assumes the form 

°^ = l-Xf3V»(x) (32) 
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with = ksT/2E m — /crjTo, the elementary magnetic flux energy E m = 

and fco = ksT* /2E m is the ratio of two characteristic energies. The dimensionless 

quantum correction parameter 



A = A 



y + & 1 



c 



An = 



hR 



h/2nCR 
kuT* 



(33) 



Remember that the Smoluchowski regime corresponds to the strong coupling limit. 
For classical systems, i.e. when the quantum correction parameter A = 0, the station- 
ary state is a Gibbs state, i.e. p(x) oc exp[— f3V (x)]. For quantum systems, due to the 
.T-dependence of the modified diffusion coefficient D\(x), the stationary state (28 1 is 
not a Gibbs state. However, it is a thermal equilibrium state. 




Fig. 1 (color online) The dimensionless current i vs the external magnetic flux x e for three values of the 
probability p of an even number of coherent electrons in the ring. Please note that one can obtain the 
persistent current for slightly asymmetric case p = 0.48 just by shifting the presented characteristic for 
p = 0.52 by x e = 0.5. It is the case of "classical" thermal fluctuations, i.e. when the quantum correction 
parameter Arj = 0. Other dimensionless parameters are: a = 0.1, Tq = 0.5, kg = 0.08. 



6 Current - flux characteristics in stationary states 

A persistent current is a periodic function of the magnetic flux with a period given 
by a single-electron flux unit </> = h/e. To what extent the current is highly sen- 
sitive to a variety of subtle effects such as an electron-electron interaction, defects, 
disorder, coupling to an environment and other degrees of freedom, it is still a topic 
of controversy and persistent discussion. Fortunately, novel techniques developed re- 
cently such as the microtorsional magnetometer |5| and scanning SQUID [6| allow 
to measure the persistent current in metal rings over a wide range of magnetic fields, 
temperatures and ring sizes. Like many mesoscopic effects, the persistent current in 
real systems depends on the particular realization of disorder and thus varies between 
nominally identical rings, cf. the term cos(nkpl) in Eq. (|9| which in practice is ran- 
dom. In Fig. 1 we depict the well-known dependence of current upon the external 
magnetic flux x e for three selected values of the probability p of an even number of 
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coherent electrons in the ring. The choice of values of p is arbitrary but shapes of the 
current are similar to those observed in experiments. 



We now focus on the influence of quantum thermal fluctuations on persistent 
currents. The deviation from the "classicality" is measured by the dimensionless pa- 
rameter A which depends on Ao (the material constant) and temperature, see Eq. ( 33 1. 
It is instructive to compare basic quantities characterizing the system. In Fig. 2 we 
show the generalized thermodynamic potential \&\{x), the modified diffusion func- 
tion D\(x) and the stationary probability density p(x) for two values of the quantum 
correction strength Ao. Three panels (a), (b) and (c) are presented for the case x e = 
(the vanishing external flux). The case Ao = corresponds to classical thermal fluc- 
tuations and Wo(x) is a bare potential V(x)/koTo. We note that the generalized ther- 
modynamic potential tf^ (x) for various A changes only slightly. On the contrary, the 
state-dependent diffusion function is a periodic function of the magnetic flux and 
possess maxima and minima. It is a radical difference to the classical case A = 
for which Do (x) = Do = fcoTo is a constant function (thin solid blue line and thin 
dashed red line in panel (c)). The maxima of D\(x) can be interpreted as a higher 
effective local temperature. They are located at x e = 1/4 mod(l/2). The impact of 
quantum corrections on the stationary probability density p(x) seems to be rather in- 
significant. One can observe a small deformation around the peak of the density: for 
lower temperature and non-zero Ao the peak becomes slightly higher and narrower 
and the tails do not diverge in the quantum case as fast as in the classical one. 



Finally, we analyze the influence of quantum thermal fluctuations on the current- 
flux characteristics i = i(x e ). It is illustrated in panel (d) of Fig. 2 and in Fig. 3. 
Two solid blue lines in panel (d) of Fig. 2 are qualitatively similar to the experimental 
curve shown in figure S6(A) in the Supporting Online Material [21 1 of the paper Q. 
We observe that in all cases of quantum thermal fluctuations the amplitude of persis- 
tent currents is reduced in comparison to classical thermal fluctuations case (Ao = 0) 
both in the paramagnetic regime (To = 0.5) and diamagnetic regime (To = 0.8). The 
parameter regime depicted in Fig. 3 is much more interesting. For T = 0.6, in the 
"classical" case, the persistent current is paramagnetic, i.e. i = r]x e with the positive 
slope r) > in the vicinity of x e = 0. It is a linear response regime where the trans- 
port coefficient (susceptibility) r) = ]imx e -^o[i{x e )/x e ]. If temperature is a Utile bit 
higher (To = 0.642), the susceptibility r/ = zero in the classical case. If quantum 
corrections are taken into account, the susceptibility 77 < , the slope of the i — x e 
curve in negative and the current becomes diamagnetic. The most interesting obser- 
vation is that the persistent current can change its character from the paramagnetic 
to diamagnetic phase and the sign of the low-field magnetic response depends on the 
level of the quantum corrections. Our detailed numerical analysis shows that the sign 
of magnetic susceptibility can easily be affected by system parameters and therefore 
is not robust against small perturbations. This is what has been observed in many 
experiments regarding the paramagnetic or/and diamagnetic persistent currents. The 
best illustration of what we state here is the response of 15 nominally identical ring 
presented in Fig. 2 in Ref. [6 |: e.g. for the ring 1 the current is paramagnetic while 
for the ring 2 it is diamagnetic. 
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Fig. 2 (color online) Four characteristics of the normal metallic ring in the absence (panels a - c) and 
presence (panel d) of the external flux are demonstrated for two values of the dimensionless temperature: 
Tq = 0.5,0.8 and for two different values of the quantum fluctuations strength: Ao = 0,0.001. The 
key for reading this plot is as the following: blue (solid) lines denote curves corresponding to the lower 
temperature To = 0.5 and the red (dashed) lines correspond to the higher temperature To = 0.8; the thin 
lines correspond to Ao = (classical thermal fluctuations) and the thick lines correspond to Ao = 0.001 
(quantum thermal fluctuations). In panel (a) we present the thermodynamic potential defined in \29\ . In 
the classical case (Ao = 0) it reduces to V(x)/koTo. One can notice only small deviations from the 
classical case when the rescaled parameter Aq is increased to 0.001. On the panel (b) we illustrate the 
corresponding stationary probability density function p(x). The bell - shaped curve in the classical limit 
is slightly deformed for non-zero Ao . The peak tends to be narrower and reach higher values and the tails 
decay slower as we increase Ao . This effect is stimulated collectively by the thermodynamic potential and 
the effective diffusion presented on panels (a) and (c), respectively. Please note that for the classical case, 
the flux dependence of D\(x) disappears and constantly equals fcoT) = 0.04 for To = 0.5 and 0.064 
for 0.8. The most significant influence of the temperature is found in panel (d), where the current-flux 
characteristics are displayed. For small external load, around flux x e = 0, the system responses in a 
completely different way for two selected temperatures. For To = 0.5 in both classical and quantum cases 
the persistent current is paramagnetic. If we, however, increase the temperature to To =0.8, the situation 
changes drastically and the susceptibility for this higher temperature becomes diamagnetic. Other rescaled 
parameters are set as the following a = 0.1, feo = 0.08, p = 0.48, e = 100. 

7 Conclusions 

In many cases and for various systems at the "intermediate" temperatures, the semi- 
classical theory is insufficient and the quantum corrections should be involved. It has 
been shown in the literature that in the strong friction limit, the quantum effects are 
restricted not only to low temperatures and therefore they should be incorporated for 
the higher temperatures as well. This is so because the quantum fluctuations, even if 
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Fig. 3 (color online) The dimensionless current i vs the external magnetic flux x e for two values of the 
dimensionless temperature To = 0.6, 0.642 and two values of the parameter Aq = (classical thermal 
fluctuations) and 0.001 (quantum thermal fluctuations). Again, the key for reading this plot is as for the 
previous figure: blue (solid) lines denote curves corresponding to the lower temperature To = 0.6 and the 
red (dashed) lines correspond to the higher temperature To = 0.642; the thin lines correspond to Ao = 
and thick lines corresponds to Ao = 0.001. The most significant influence of the 'quantum parameter' 
Ao is found for the rescaled temperature To = 0.6, where the current i changes its behavior from para- 
magnetic to diamagnetic one just by adjusting Ao from to 0.001. Moreover, for the presented set of 
the system parameters we can find characteristic cross-temperature at To = 0.642 where the magnetic 
susceptibility is zero in the classical case and diamagnetic in quantum, see red (dashed) line for details. 
Rescaled parameters are set as the following a = 0.1, fco = 0.08, p = 0.48, e = 100. 

reduced for one variable, are enlarged for the conjugate variable. The dynamics as 
well as the stationary states in this regime can be modeled by the quantum Smolu- 
chowski equation. In other words, the quantum non-Markovian stochastic process is 
approximated by the classical Markovian process with the modified, state-dependent 
diffusion function. 

The role of the quantum corrections on the current-flux characteristics is ad- 
dressed in this work. A general conclusion is that the quantum thermal fluctuations re- 
duce the amplitude of the persistent currents: the current amplitude is always smaller 
than the corresponding "classical" one. In the quantum case, the diffusion constant 
becomes a periodic function of the magnetic flux. Maxima and minima of the diffu- 
sion function can be interpreted in terms of the higher and lower local temperature. 
There are parameters regimes where the system response changes the character from 
the paramagnetic to diamagnetic, when the quantum effects of thermal fluctuations 
increases. It would be interesting to extend the current study by including the time- 
dependent drivings modeled by the time-periodic magnetic fields. One could expect 
novel transport phenomena like a negative susceptibility which for Brownian parti- 
cles corresponds to negative mobility ll22l or negative conductances 1 14 1. 
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